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Abstract. Linear polarisation in gamma-ray burst prompt emission is an important diagnostic with the potential to significantly 
constrain models. The spectrometer aboard INTEGRAL, SPI, has the capability to detect the signature of polarised emission 
from a bright y-ray source. GRB 041219a is the most intense burst localised by INTEGRAL with a fluence of 5.7 xlO"'* 
ergs cm"- over the energy range 20 keV-8 MeV and is an ideal candidate for such a study. Polarisation can be measured using 
multiple events scattered into adjacent detectors because the Compton scatter angle depends on the polarisation of the incoming 
photon. A search for linear polarisation in the most intense pulse of duration 66 seconds and in the brightest 12 seconds of 
GRB 041219a was performed in the 100-350 keV, 100-500 keV and 100 keV-1 MeV energy ranges. It was possible to divide 
the events into six directions in the energy ranges of 100-350 keV and 100-500 keV using the kinematics of the Compton scatter 
interactions. The multiple event data from the spectrometer was analysed and compared with the predicted instrument response 
obtained from Monte-Carlo simulations using the GEANT 4 INTEGRAL mass model. The;^'^ distribution between the real and 
simulated data as a function of the percentage polarisation and polarisation angle was calculated for all three energy ranges. The 
degree and angle of polarisation were obtained from the best-fit value of x~- A weak signal consistent with polarisation was 
found throughout the analyses. The degree of linear polarisation in the brightest pulse of duration 66 s was found to be 63^3^% 
at an angle of lO^W degrees in the 100-350 keV energy range. The degree of polarisation was also constrained in the brightest 
12 s of the GRB and a polarisation fraction of 96;^;^^% at an angle of 60;^|^ degrees was determined over the same energy range. 
However, despite extensive analysis and simulations, a systematic effect that could mimic the weak polarisation signal could not 
be definitively excluded. Our results over several energy ranges and time intervals are consistent with a polarisation signal of 
about 60% but at a low level of significance (~ 2cr). The polarisation results are compared with predictions from the synchrotron 
and Compton drag processes. The spectrum of this GRB can also be well fit by a combined black body and power law model 
which could arise from a combination of the Compton and synchrotron processes, with different degrees of polarisation. We 
therefore conclude that the procedure described here demonstrates the effectiveness of using SPI as a polarimeter, and is a 
viable method of measuring polarisation levels in intense gamma-ray bursts. 
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1. Introduction 

Polarisation is a powerful tool for investigating emission pro- 
cesses in long gamma-ray bursts (GRBs). The link between the 
y-ray production mechanism and the degree of linear polarisa- 
tion can be exploited to constrain models. 

Long gamma ray bursts are linked to the collapse of a mas- 
sive star which forms a rapidly rotating black hole. For a recent 
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review of GRBs, see Meszaros (l2006l) . In addition, a large or- 
dered magnetic field m ay be induced by the angular momentum 
of the accretion disk dZhang & MeszarosI l2004t iPiranI 120041) . 
Energetic outflows develop which are beamed perpendicular to 
the accretion disk and along the black hole's rotation axis. An 
observer close to the jet axis will detect a GRB. Polarisation 
is generally associated with an asymmetry in the way that the 
material is viewed. The asymmetry can be attributed to a pref- 
erential orientation of the magnetic field or to inverse Compton 
scattering. The polarisation mechanisms are discussed in more 
detail in ^ 



The reported detection of significant polarisation 
(Oj = 80 + 20% in the energy range 15-2000 keV) in 
GRB 021206 dCoburn & Boggs, ,2003,) using the RHESSI 
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spacecra ft led to many publications examining the results 



(Rutled ge&Foxl l2004t IWigger et alj 12004 iBoggs & Coburn 



[20031) and the mechanisms for producing l a rge polarisation 



(e.g. IShaviv & Daii Il995l: iNakar et all 120031: IWaxmalJ 12003 



20071). The 



GranotI 120031: iLazzati et all 12004 iDado et al 
RHESSI results highlighted the importance of correctly 
evaluating the systematic effects, which m ay mimic a polar- 



isation signature. A recent novel attempt dWillis et alj l2005h 



involved analysing the Earth's albedo flux seen by BATSE 
for GRB 930131 and GRB 960924, where the lower limits 
of polarisation were found to be 11 , > 35% and 11 , > 50% 
respectively. These figures can only be considered as lower 
limits due to systematic effects, including natural anisotropics 
in the Earth's albedo flux and possible limitations in the 
GEANT 4 code at the time the simulation was run. 

The dominant mode of interaction for photons in the energy 
range of a few hundred keV is Compton scattering. Linearly 
polarised y-rays preferentially scatter perpendicular to the in- 
cident polarisation vector, resulting in an azimuthal scatter an- 
gle distribution (ASAD) which is modulated relative to the dis- 
tribution for unpolarised photons. The sensitivity of an instru- 
ment to polarisation is determined by its effective area to scat- 
ter events and the average value of the polarimetric modulation 
factor, Q, which is the maximum variatio n in azim uthal scat- 
tering probability for polarised photons (iLei et al ] 11997) . The 
value of Q is given by 



dcrj_ — dcr\\ 
dcr± + do-\\ 



(1) 



where dcrj^, dcry are the Klein-Nishina differential cross- 
sections for Compton scattering perpendicular and parallel to 
the polarisation direction, respectively. Q is a function of in- 
cident photon energy, E, and the Compton scatter angle, 0, be- 
tween the incident and scattered photon directions. For a source 
of count rate S and fractional polarisation Ilj, the expected 
ASAD is given by: 



^=(|;:)[l-Gn.cos2(0-77)] 



s 



(2) 



where (p is the s cattering angle, and rj is the polarisation angle 
dLei et al.lll99 7'). This equation yields a 180° modulated curve 
when fit to polarised data, where rj represents the minimum 
angle of the modulated distribution and gives the direction of 
the polarisation vector. 

1.1. SPI as a polarimeter 

SPl is not optimised to act as a polarimeter, but because of its 
detector layout, geometry and thick detector plane, the mod- 
ulation from a polarised flux can be measured through multi- 
ple scatter events in its detectors. Kalemci et al.. (2004) found 
that it is possible to measure polarisation in a moderately 
bright GRB in the field of view of SPI if the GRB is on- 
axis. GRB 041219a had a fluence of 5.7 xlO""* ergscm"^ and 
a peak flux of 1.84 xlO"^ ergs cm"^ ^-i (20 keV-8 MeV) 
at an off-axis angle of 3.2° and is the most intense burst de- 
tected by INTEGRAL, so would appear to be an ideal candidate 



( McBreen et alj2006l ). Detailed Monte-Carlo simulations built 
with the GEANT 4 toolkit can be used to predict the response 
of SPI to a polarised flux. A comparison of the data and simu- 
lations enables a determination of the polarisation strength and 
angle. 




Fig. 1. The numbering system used for the Germanium SPI de- 
tectors. The 6 directions used in the polarisation analysis are 
shown along with the y- and z-coordinate axes of the space- 
craft. The X-axis is normal to the detector plane. 



2. The Spectrometer aboard INTEGRAL 

The European Space Agency's International Gamma-Ray 
Astrophysics Laboratory, INTEGR AL, was launched on 17 
October 2002 dWinkler etaDl2003 f). It consists of two coded 
mask y-ray instruments, the spectrometer (SPI) and the imager 
(IBIS). The instruments are coaligned so that data is taken by 
all instruments in one pointing. 

SPI consists of 1 9 hexagonal germanium (Ge) detectors 
dVedrenne et al ]|2003h . arranged to minimise the volume of the 
array and the space between each detector (Fig. [T]!. The de- 
tectors cover the energy range 20 keV-8 MeV with an energy 
resolution of 2.5 keV at 1.3 MeV. Each detector is 6.9 cm in 
height, with a centre to centre distance of ~ 6 cm between adja- 
cent crystals. A coded mask is located 1.71 m above the detec- 
tor plane for imaging purposes, giving a 16° corner-to-corner 
field of view. The sensitivity of SPI (~ 5 xlO"^ photons cm"^ 
s"^ keV"') is limited by the instrumental background, which 
consists mainly of cosmic rays impinging on the det ectors and 
the secondary particles created by their interaction ('jean et al] 
,200 3: Wei denspointner et al. 2003 ). The background can be de- 
termined by averaging the count rate over a long period of time 
during the science window, and subtracting this average from 
the raw count rate. The background is significantly reduced by 
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the presence of an anti-coincidence shield made from BGO 
crystals surrounding the Ge detectors. 

The operating mode of SPI is based on the detection of 
events from the Ge detectors which are not accompanied by 
a corresponding detection in the anti-coincidence shield. The 
events are separated into single events (SE) where a photon de- 
posits energy in one detector, and multiple events (ME) where 
the photon deposits energy in two or more detectors. All events 
are processed by the Digital Front End Electronics (DFEE), 
which provides event timing and classification. SPI operates in 
photon-by-photon mode, which produces photon packets (80 
packets/8 s) containing all of the non-vetoed events and scien- 
tific housekeeping packets (5 packets/8 s) including the event 
counters which are used to generate lightcurves. 

Detectors 2 and 17 ceased to function on December 6, 
2003, and July 17, 2004 respectively. The failure of these detec- 
tors results in a decrease of the effective area of the instrument 
to about 90% of the original area for SEs. It is reduced to ~ 
75% for MEs, because the number of pseudo detectors (i.e. the 
adjacent detector pairs used to measure multiple events) drops 
from 84 to 64. 

3. Model simulation 

The advent of fast computing clusters has made difficult 
computational tasks such as the prediction of instrument re- 
sponse to polarised flux more feasible. A computer model of 
the INTEGRAL spacec raft written in the GEANT 4 toolkit 
(lAgostinelli et al ] l2003h was used for simulations in this work. 
This model was develope d from the GEAN T 3 INTEGRAL 
Mass-Model (TIMM) CFerguson et al.1 l2003b originally used 
to assess the background recorded by the instruments onboard 
INTEGRAL. The model contains an accurate representation of 
the SPI instrument, including the mask and veto elements. The 
rest of the spacecraft is modelled to a much lower level of de- 
tail. Average densities and simplified geometries are used for 
areas of the spacecraft positioned at larger distances from the 
detectors since they will not have a large effect on the on-axis 
gamma-rays. 

The GEANT 4 toolkit contains all the physics necessary 
to allow the tracking of photons and particles through a mod- 
elled geometry. The software consists of a series of random 
number generators to calculate the probability of an interaction 
occurring in a material. As with any program, the simulation 
is dependent on the coding of the interaction in the software. 
Mizuno et al.l (120051) reported an incorrect implementation of 
the polarised Compton and Rayleigh scattering processes in the 
GEANT 4 code. This error caused the azimuthal modulation 
due to the polarisation to be lower than expected. When the 
appropriate correction was applied, an increase of ~15% was 
seen, even in the higher energy regime used in our simulations 

(Fig.ia. 



3.1. Simulating GRB 041219a 

Gamma-ray photons were directed into the model geometry 
from a plane surface in the direction of the GRB, 3.08° from 
the INTEGRAL x-axis and 63.95° from the INTEGRAL z-axis. 



Comparison Dt modjlalion between original and modiried GEANT4 code 
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Fig. 2. Simulated modulation due to Compton Scattering in 
a test geometry. The solid line gives the original GEANT 4 
Compton scatter code and the dashe d line gives the Comp ton 
scatter code using modifications from lMizuno et alj (l2005h . 




Fig. 3. The mask elements (yellow) overlaying the detectors 
(blue), as viewed from the direction of the incoming GRB pho- 
tons generated from the simulations. Fig. [T| shows the number 
allocated to each detector. 



simulating the incoming ffux from a source at infinity from the 
same d irection relative to the spacecraft as the GRB. The Band 
model jBand et al. I ll993h parameters for the main peak of the 

1 +0.08 
-0.21' 



burst of duration 66 seconds (a - -1.50 ^o o6' ■ 

keV) were used to create the spectrum (see ^4.2i . 



f310 1 



Eo = 568 : 

This 66 s interval was selected to maximise the source counts. 

For each simulation run, the polarisation angle of the pho- 
tons was set between 0° and 180° in 10 degree steps, and the 
polarisation fraction was set to 100%. There was also one run 
for a beam of unpolarised photons. Only polarisation angles 
between 0° and 180° were simulated due to the symmetry of 
the system and the difficulty in separating the scattering direc- 
tions between the pixels. The effect of the spectral shape on 
the level of polarisation was simulated, and it was found that 
the simulated polarisation signal depended very weakly on the 
spectral parameters. The secondary photons produced by the 
multiple events scatter more in the forward direction at higher 
energies, causing the azimuthal modulation to drop sUghtly. At 
lower energies the multiple events are less likely to occur due 
to the photoelectric effect dominating the scatter processes. 
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The simulations produced a list of all the interactions that 
occurred in the sensitive volumes of the model (Ge detectors 
and BGO shield). These data were then converted into an event 
list, for comparison to the real SPI data. Initially the interac- 
tions were summed, so that the energy deposits correspond to 
the total energy deposited for an event in each of the sensitive 
volumes. These deposits were then filtered according to the en- 
ergy thresholds of the detectors (~20 keV) and veto (~80 keV). 
After subtracting the vetoed events, the event list was separated 
into single events (where the photon is detected in one pixel) 
and multiple events (where the photon is detected in multiple 
pixels). This process produced the final list of events to analyse 
and compare to the real data. The unpolarised simulation data 
was combined with the polarised simulation data, allowing the 
percentage of polarisation to be changed for each angle. 

4. SPI Data Analysis 

GRB 041219a was det ected by IBAS a t 01:42:18 UTC on 
December 19* 2004 (iGotz et al.1 |2004 at a location of 
00h24m25.8s, +62° 50' 05.6" close to the axis of the detec- 
tor. 

4.1. Temporal Analysis 



3500 
3000 
2500 
2000 
1500 
1000 
500 


-500 




100 200 300 400 500 

Time from 40s before to 520s after 01 :42:1 8 UTC 

Fig. 4. Background-subtracted single event lightcurve of 
GRB 041219a, summed over all SPI detectors in the energy 
range 20 keV-8 MeV. The vertical solid lines mark the start 
and end of the 66 second emission phase (To = 261 s to To = 
327 s). The vertical dashed lines mark the start and end of the 
brightest 12 seconds of the burst (To = 276 s to To = 288 s). To 
is the IBAS trigger time (01:42:18UTC). 



GRB 041219a consisted of an initial precursor-type pulse, 
followed by a quiescent period lasting approximately 200 s, 
before the main emission beginning at ~ 250 s post-trigger. 
An image of the coded mask as seen from the direction of 
the incoming GRB photons was obtained from the simulations 
(Fig. Is). The background-subtracted single event lightcurve 
summed over all of the detectors was generated and is shown 
in Fig. m The mask almost completely obscured three of the 
detectors (12, 3, 0), and partially obscured five more (4-6, 8, 
13) (Fig.|6]l. Also, detectors 2 and 17 are no longer functioning 



(and were not included in the analysis). However, the GRB can 
be clearly recognised in at least nine of the SB lightcurves in 
Fig.|6] 



1000 

channel energy (keV) 




Energy (keV) 



Fig. 5. (a) Band model fit to GRB041219a and residuals 
for the 66 s emission phase (Fig. |4|i. (b) vFy spectrum of 
GRB 041219a. The Band model parameters are o- = -1.50, 
yS= -1.95 and £o = 568 keV. 



4.2. Spectral Analysis 

The spectrum of GRB 041219a was extracted us ing specific 



GRB tools from the Online So ftware Analysis (IDiehl et al 



2003; Skinner & Connell 2003) version 5.0 available from 



the INTEGRAL Science Data Centre. GRB 041219a is the 
brightest burst localised by INTEGRAL with a peak flux of 
43phcm"^s"' (20keV-8 MeV). The spectrum of t he burst and 



sub-in tervals were well fit by the Band model (iBand et al 



19931), although the parameters of the spectrum evolved during 
the burst. A detailed discussion of t he spectral and tempor al be- 
haviour of this burst is available in McBreen et al. I ( I2OO6I) . The 
most intense emission pulse of duration 66 s (indicated by the 
solid lines in Fig.|4]i was selected for polarisation analysis. The 
photon indices, a and yS, for the emission phase used to calcu- 
late the polarisation were -1.50 ^^jj ^g and -1.95 ;^q 2[ respec- 
tively. The break energy Eo was 568 ^j^j keV. The spectra are 
shown in Fig.|5] The peak energy, iSpeak, is given by (2 + a)y. Eq 
and the value of fipeak in the interval of the main emission phase 
is 284 keV. In addition, the polarisation analysis was per- 
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Fig. 6. The layout of the 19 detectors of SPI with single event lightcurves of GRB 041219a showing the variation in count rate 
per detector. The horizontal and vertical axes give the time and count rate in each detector respectively. The detector number 
is indicated in the corner of each lightcurve. Detectors 2 and 17 are not in operation. The detectors with high count rate were 
unobscured or partially obscured by the mask (Fig.|3]l. 



formed for the brightest 12 s of the 66 s interval to determine 
the polarisation over the duration of this intense pulse. 

It is interesting to note that the spectrum of GRB 041219a 
was equally well fit by a combina t ion of a black body and 



power law model (iMcBreen et al. I2OO6I) . iFan et al.1 tOOSh 



also found that the early optical and infrared emission from 
GRB 04 1 2 1 9a can be modelled as the superposition of a reverse 
and a forward shock component. The ejecta are magnetised to 
a small extent, which may be d ue to magnetic fi eld generation 
during the internal shock phase. iFan et al. I (l2005h predicted that 
the internal shock emission was very likely to be linearly po- 
larised. 



5. Polarisation Analysis 

There is no positional resolution within the SPI detectors and 
so it is not possible to determine the exact position of the in- 
teraction within an individual detector. Centre-to-centre inter- 
actions are assumed for multiple events. According to the sim- 
ulations, this will introduce an uncertainty on each angle of 
~29°. Below 511 keV, the incoming photons predominantly 
Compton scatter from the detector with the lower energy de- 
posit to the higher one (IKalemci et al ]|2004 . Thus 6 directions 
of scatter can be distinguished. For higher energies, the order 
of energy deposition does not distinguish between anti-parallel 
directions and the number of directions is limited to 3. To en- 
able a larger energy range from 100 keV-1 MeV to be investi- 
gated, the analysis was also performed in 3 directions. 



5.1. Method 

The analysis procedure was carried out, starting with the raw 
data from SPI, as follows: 



1. All interactions between detectors (double events) during 
the defined time intervals were selected. 

2. Only double events that occurred in adjacent detectors were 
accepted. 

3. The list of events was calibrated (using the spi_gain_corr 
tool from OSA 5.0) to convert the original channel number 
to energy in keV. 

4. All interactions with less than 30 keV deposited per detec- 
tor were rejected. Relatively few photons above 100 keV 
will lose < 30 keV in Compton scatter interactions. 

5. Coincident pairs whose combined energies lie in the 100- 
350 keV range were selected. Up to 511 keV, incoming 
photons predominantly scatter from the detector with the 
lower energy deposit to the higher one allowing 6 separate 
directions to be determined. 

6. The scatter pairs were divided into 6 different directions 
(0-300 degrees) and the total number of events in each di- 
rection were determined (Fig.[T]i. 

7. Background events (using the same selection process) were 
selected from intervals in the same science window before 
the GRB occurred (Table [TJ, since the emission continued 
up to the end of the science window. The scaled background 
was then subtracted from the multiple event list. The set of 
19 detectors was also divided into 4 quadrants and the total 
background was calculated separately for each quadrant to 
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Table 1. Time intervals used for the background determination 
in 6 directions. 



Observation 


Start Time (UTC) 


End Time (UTC) 


Background 1 


01:30:00 


01:31:06 


Background 2 


01:35:00 


01:36:06 


Background 3 


01:36:30 


01:37:36 


Background 4 


01:38:00 


01:39:06 


Total (sec) 




264 



ensure that there were no biases in any specific direction or 
systematic eff'ects. 

The analysis was carried out for 6 directions in the energy 
ranges 100-350 keV and 100-500 keV and over two sepa- 
rate time intervals (Fig. |4| as described above. The analysis 
was also performed for 3 directions in the 100-350 keV, 100- 
500 keV and 100 keV-1 MeV energy ranges to compare the 
values obtained from both methods. The number of multiple 
events between 100-350 keV, 100-500 keV and 100 keV- 
1 MeV were 860, 1218 and 1876 respectively for the 66 s time 
interval, and the total number of simulated events was ~ 10^ 
per energy range. The simulated and real data sets were scaled 
by the total number for all directions to ensure that the compar- 
isons between both types of data were valid and anisotropics in 
the response due to the mask and the two inoperative detectors 
were taken into account. 

Each second of data during the 66 seconds was also anal- 
ysed separately using the standard OSA software. It was ob- 
served that approximately 30 seconds into the brightest portion 
of the burst, the live time per second of each detector dropped 
dramatically to about half of its original value due to the high 
data rate and telemetry limitations (Fig.|7]). The result was that 
almost half of the multiple events during this period were lost, 
and so it was necessary to reduce the ME background to take 
this loss into account (Fig.[8]l. The analysis was carried out for 
the brightest 12 seconds of the pulse (To = 276 s to To = 288 s) 
before the event rates were significantly afifected by packet loss 
to check the method and to avail of a higher signal to noise. The 
two sets of results for the 12 second and the 66 second intervals 
could then be analysed separately and compared. 

The multiple event rate between detectors on opposite sides 
of the array was examined to determine the random rate be- 
tween non-adjacent detectors and to investigate if the GRB 
had sub-microsecond variability (i.e. if events were deposited 
in a shorter interval than the 350 ns coincidence time window). 
It was observed that even between detectors with high single 
event count rates (e.g. detectors 10 and 15) in the 66 second 
interval, the average multiple event rate was approximately 1 
count over the 66 second duration. This result agrees with the 
expected random rate and excludes sub-microsecond variabil- 
ity in GRB 041219a. 

The event distri bution is highly dependent on SPI's geome- 
try jLei et al Since there are inhomogeneities in the de- 
tector layout (e.g. inoperative detectors and detectors covered 
by the coded mask), the Q distribution will also be distorted. 
This distribution as a function of polarisation angle was simu- 
lated and taken into account when estimating an average value 



0.9 
0.8 
0.7 

■I 0.6 

> 
_j 

0.5 



0.4 



0.3 




260 



270 



hmmk 

280 290 300 310 
Time since trigger (s) 



320 



330 



Fig. 7. The average live time for each SPI detector per second 
over the most intense phase of emission of the GRB, showing 
a marked decline ~30 seconds into the pulse. 





90 - 




80 - 




70 - 




60 - 


c 


50 - 


ZJ 
O 




o 


40 - 




30 - 




20 '.f 




10 - 







Corrected MEs _ 






Original MEs 









270 



280 290 300 
Time since trigger (s) 



310 



320 



Fig. 8. The original multiple event lightcurve (dashed) for all 
operative detectors for the 66 second interval, and the multiple 
event lightcurve (solid line) after the dead time correction was 
made. The vertical dashed lines indicate the 12 seconds used in 
the analysis. 



of Q. From our simulations, we estimated the average modula- 
tion factor Q for 100% polarisation to be 24 ± 7%, in agreement 
with the calculations of iKalemci et al.l (120041) . 



6. Results 

The 100% polarised and 0% polarised data obtained from the 
Monte-Carlo simulations for each scatter angle were combined 
to create a partially polarised signal with varying degrees of po- 
larisation. The fitting routine compared the real data with the 
partially polarised simulated data. The percentage polarisation 
was varied from 0% to 100% in steps of 10% and the angle was 
varied from 0° to 180° in 10° intervals. The real data were com- 
pared with the simulated data and the value ofx^ calculated for 
a range of angles and percentages of polarisation. These values 
were used to generate significance level contour plots (Figs.|9] 
and [Tol l, which gave a minimum at the angle and percentage 
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Fig. 9. Contour plots of the percentage polarisation as a function of the polarisation angle for the six scatter directions (0° - 360°) 
in the 12 second interval, showing the 68%, 95% and 99.7% probability contours. The plots indicate a non-zero value for the 
level of polarisation, (a) refers to the energy range 100-350 keV and (b) refers to the energy range 100-500 keV. 



of polarisation that most closely matched the real data. The re- 
sults of the fitting procedures are given in Table |2l which lists 
the percentage polarisation and the angle for the 12 second and 
66 second time intervals in the energy ranges 100-350 keV, 
100-500 keV and 100 keV-1 MeV. The eiTors quoted for the 
percentage and angle of polarisation are 1 cr for 2 parameters of 
interest. 

Fig. |9] shows the contour plots obtained by comparing the 
real and simulated data for the six scatter directions in the 12 
second interval. Fig. [TOl shows the corresponding contour plots 
for the three scatter directions in the 12 s and 66 s intervals. 
The contour plots for the 66 second interval for the six scatter 
directions are not shown, because the best fit probability indi- 
cates that the model was not a good fit to the real data, and a 
68% probability contour could not be generated. The contour 
plots indicate a non-zero value for the level of polarisation in 
all of the time intervals and energy ranges studied. 

Eight of the ten cases listed in Table|2]indicate that the per- 
centage of polarisation is greater than 50%. The best fit prob- 
ability that the simulated values match the real data is greater 
than 99.8% in the 12 second interval for the three scatter di- 
rections in the 100-350 keV energy range, corresponding to 
a percentage polarisation of 96^4^% at an angle of 60^[4 de- 
grees (Fig. [To] (a)). The best fit probability in the 66 second 
interval is greater than 98% for the three scatter directions in 
the same energy range, corresponding to a percentage polarisa- 
tion of 63+^^% at an angle of 70+1^ degrees (Fig.[l0](b)). The 
polarisation angles are consistent in all cases with a value be- 
tween 60° and 70°. Despite extensive analysis and simulations, 
we could not exclude a systematic effect that could mimic the 
weak polarisation signal. In some cases, the percentage polar- 
isation exceeds 100% when the errors are taken into account. 
This is due to the poor signal to noise of the data and possible 
systematic instrumental effects. 

The weighted mean level of polarisation was calculated for 
each time interval separately from the values listed in Table |2] 
where the best fit probability that the model was a good fit to 
the real data was greater than 90%. The level of polarisation 



for the 12 second time interval was 76 + 40% at an angle of 
67^}^ degrees and the level of polarisation for the 66 second 
time interval was 43 + 25% at an angle of 70;^[g degrees. The 
weighted mean for all cases listed in Table|2]was determined to 
be 60 + 35% for an angle of 68 + 15 degrees. 



7. Discussion 

The results obtained from our simulations and analysis are con- 
sistent with linear polarisation at about the 60% level (~ 2cr) 
at an angle of ~ 70°. It is possible that the percentage polari- 
sation varies with energy, angle and time over the duration of 
the burst. However, the levels of polarisation measured during 
the brightest 12 seconds of GRB 041219a and the brightest 
66 second pulse are consistent at the ~ 2cr level, indicating 
that there is no major variation in polarisation during the in- 
tense 66 second pulse. It is unlikely that a burst brighter than 
GRB 041219a will be detected by INTEGRAL. A GRB of sim- 
ilar fluence but over a shorter time interval may produce better 
statistics. Another po ssibility is a spectrall y harder burst, sim- 
ilar to GRB 941017 dGonzalez et al.ll2003h . which would pro- 
duce more multiple events in the MeV energy range and thus 
cre ate a strong polarisati on signature. 



Kalemci et al.l (120061) have independently analysed the SPI 



data for GRB 041219a, with simulations performed using the 
MGEANT code rather than the GEANT 4 code used here. By 
fitting the azimuthal scatter angle distribution of the observed 
data over the 6 directions, we obtain results consistent with 



Kalemci et al.l (12006 ) in both magnitude and direction, within 



the Umits given by the large error bars. However, the more 
complete analysis presented here compares the observed data 
to various combinations of the simulated polarised and unpo- 
larised dat a (Figs. |9] and [lOl Tab le |2]|. We agree with the con- 
clusions of Kalemci et al. |(|2006') that there is a possibility that 
instrumental systematics may dominate the measured effect. 

There are a number of different methods of measuring 
polarisation using the INTEGRAL instruments. For example, 
iMarcinkowski et al.i (,2006.) described a new method of using 
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Fig. 10. Contour plots of the percentage polarisation as a function of the polarisation angle for the three scatter directions (0°, 
60° and 120°), showing the 68%, 95% and 99.7% probability contours, (a) and (b) refer to the energy range 100-350 keV, (c) 
and (d) to the energy range 100-500 keV, and (e) and (f) to the energy range 100 keV-1 MeV. The plots on the left represent the 
12 s interval, and the plots on the right represent the 66 s interval. 



the IBIS instrument in Compton mode to detect and analyse 
an intense burst that was outside the coded and partially coded 
field of view of IBIS. GRB 030406 was well detected through 
the shield using this method. Since IBIS consists of two lay- 
ers of detector arrays ( Ubertini et al .1120031) . Compton scatter- 
ing can be used to detect the events which interact in one 
layer and scatter into the second layer. The Compton mode 
determines the energy deposit and position of the event in 



each array. Therefore, it may be possible to extend this tech- 
nique to measure the polarisation fraction of a spectral ly hard 
GRB as well as the spectral and temporal parameters. iFinger 
(2006) is also investigating the possibility of using th e IBIS 
Compton mode to search for GRB polarisation. iHaidas (2006) 
is using the RHESSI spectrometer to set instrumental limits on 
the minimum detectable polarisation for several sources, in- 
cluding GRBs. It should be noted that the BAT detector on 



McGlynn et al.: Polarisation Measurements of GRB 041219a 



9 



Table 2. Table of results from;\f^ fitting of real and simulated data. The columns from left to right list the duration of the interval, 
the polarisation percentage, angle and best-fit probability that the model simulations matched up with the real data, the energy 
ranges analysed over six directions (columns 3 and 4) and the energy ranges analysed over three directions (columns 5-7). The 
errors quoted are Icr for 2 parameters of interest. 





Polarisation 


6 Directions (Fig.[9]( 
100-350 keV 100-500 keV 


3 Directions (Fig.llOll 
100-350 keV 100-500 keV 100 keV-1 MeV 


12 second 
interval 


Percentage (%) 
Angle (°) 
Probability (%) 


98 ± 53 71 +^2 
60 70 *_fi 
87.2 93.4 


96 70 ± 37 68 ± 29 
60+'^ 7Q+15 TO^'"* 

""-14 '"-14 '"-10 

99.8 99.5 95.9 


66 second 
interval 


Percentage (%) 
Angle (°) 
Probability (%) 


70 ±20 52 ±11 
70 70 ± 5 
18.4 36.9 


63 49 ±24 26 ± 20 
1Q*\\ 

98.0 95.5 97.3 



SWIFT dGehrels et al.ll2004 is not configured for polarisation 
measurements of GRBs. However, a number of missions have 
been pro posed specifically t o measure GRB polarisation e.g. 
POLAR dProduit et al] |20Q5') and XPOL feosta et alj|2006l) . 



The spectra of GRB 041219a have been well fit by both 
the Band m odel and a combinatio n of a bla ck body plus power 
law model jMcBreen et al.ll2006i) . Recently lRvd3 (l2005l) stud- 
ied the prompt emission from 25 bright GRBs and found that 
the time resolved spectra could be equally well fit by the 
black body plus power l aw model and with the Band model. 
Rees & Meszaro sl (l2005h suggested that the Epeak in the y-ray 



spectrum is due to a Comptonised thermal component from the 
photosphere, where the comoving optical depth falls to unity. 
The thermal emission from a laminar jet when viewed head- 
on would give rise to a thermal spectrum peaking in the X-ray 
or y-ray band. The resulting spectrum would be the superpo- 
sition of the Comptonised thermal component and the power 
law from synchrotron emission. Unfortunately, the polarisation 
measurements of GRB 041219a are not sensitive enough to de- 
tect the change in polarisation that might result from the com- 
bination of the Compton and synchrotron processes. 

A significant level of polarisation can be produced in GRBs 
by either synchrotron emission or by inverse Compton scatter- 
ing. The fractional polarisation produced by synchrotron emis- 
sion in a perfectly aligned magnetic field can be as high as 

- (p + l)/{p + 7 /3) where p is the power law index of 
the electron distribution. Typical values of p - 2-3 correspond 
to a polarisation of 70-75%. An ordered magnetic field of this 
type would not be prod uced in s hocks but could be advected 
fro m the cen tral engine jGranot & Konigl 2003; Granot 2003i 
Lvutikov et al.ll2003l) . 

Another asymmetry capable of producing polarisation, 
comparable to an ordered magnetic field, involves a jet with a 
small opening angle that is viewed slightly off-axis dWaxman 
I2OO3V A range of magnetic field configurations have been con- 
sidered (iSari,,1999.: ,Ghisellini & Lazzati..l999: .Granot ,2003.: 



Nakar et alJl2003t iFan et alJl2005h . The intensity distribution 
and maximum polarisation of the jet are modified if the pitch 
angle distribution of the electrons is not isotropic , but biased 
towards the orthogonal direction dLazzati I2OO6'). The more 
anisotropic distribution produces larger net polarisation. For 
broader jets, only a small fraction of random observers would 
detect a high level of polarisation. 

Shaviv & DaJ d 1995b and lOar & de Rululal d2004 have 



pointed out that polarisation is a characteristic signature of the 
inverse Compton process. This mechanism was also considered 
in the framework of an ensheathed fireball tEichler & Levinson 



2003h . Compton Drag (CD) emission is produced when an 



ionised plasma moves relativistically through a photon field. 
A fraction of the photons undergo inverse Compton scattering 
on relativistic electrons and have their energies increased by 
~4y^ where y is the electron Lorentz factor, and under certain 
circumstances the scattered photons have high polarisation. 

Lazzati et alJ d20()4|) considered CD from a fireball with an 



opening angle comparable to the relativistic beaming. The po- 
larisation is lower than that from a point source because the 
observed radiation comes from different angles. In the fireball 
model, the fractional polarisation emitted by each element re- 
mains the same, but the direction of the polarisation vector of 
the radiation emitted by different elements within the shell is 
rotated by different amounts. Th is can lead to effectiv e depo- 
larisation of the total emission (^ Lvutikov et al. 2003), which 
is not observed in GRB 041219a. A lower level of polarisa- 



tion h as recently been predicted for X-ray flashes dDado et al 



2007). 



Lazzati et alJ d2004ft calculated the polarisation as a func- 



tion of the observer angle for several jet geometries, and 
showed that a high level of polarisation can be produced if the 
condition FOj < 5 is satisfied, where F is the Lorentz factor 
of the jet and Oj is the opening angle of the jet. In the case of 
GRB 041219a, it is possible to estimate the values of F and 9j 
in the following way. GRB 041219a is estimated to have a red- 
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shift of z ~ 0.7 using the Yonetoku relationship (lYonetoku et al 

ergs 



2004). The fluence from 20 keV to 8 MeV is 5.7 x lO""* 



cm ^, yielding a value of ~1 0^"^ ergs for the total isotropic emis- 
sion dMcBreen et al. I l2006l ). The standard beaming co rrected 
energy for GRBs is £ = 5 x 10^° ergs dFrail et all I2OOI). 
Combining this information with the total isotropic emission 
yields a value of 0j ~ 2.5° (0.044 rad). The Lorentz factor of 
the fireball can be obtained from the redshift corrected peak 
energy Ep.ak (Ep^ak = 483 keV for GRB 041219a) by the rela- 
tionship 

Ep,,, ^lOrHT (3) 

where T~10^ K is th e black body spectrum of the photon field 
( Lazzati et alJ 2004 ). The computed value is F ~ 75, yielding 
the result: 

TOj ~ 3.3 (4) 

The small value of FOj shows that it is possible to have 
polarisation of ~6Q% in GRB 041219a and also produce the 
lower limit to the value s of the polarisation for two BATSE 
GRBs (I Wilhs et alJ2005h and the value of 41+ ^^ % obtained by 
RHESSI for GRB 021206 dWigger et alJl2004 . 

Synchrotron radiation from an ordered magnetic field ad- 
vected from the central engine and Compton Drag are both 
good explanations for a significant level of polarisation. It 
should be possible to distinguish between the two emission 
mechanisms. Only a small fraction of GRBs should be highly 
polarised from Compton Drag because they have narrower jets, 
whereas the synchrotron radiation from an ordered magnetic 
field should be a general feature of all GRBs. Another possible 
distinction between the two processes involves the optical flash 
because the Compton Drag radiation should be less polarised 
than synchrotron radiation. 

A small but significant degree of linear polarisation 
was discovered in the optical afterg low of GRB 990510 
dCovino et al. I ll 9991: IWiiers et al Since then, there have 

be en a number of other detections of polaris ation in afterglows 
(eglCovinoetal .12004 iGorosabel et al E004t) . The polarisation 
is observed to be at about the 1-3% level and is reasonably 
constant when asso ciated with a smooth afterglow lightcurve 
dCovino et al ]l2003h . The polarisation can vary in direction and 
degree on a time scale of hours if there are de viations from the 
smooth power law decay dGreiner et al.ll2003h . For a review of 
the levels of asymmetry needed to provide a polarisation signal 
in the prompt and afterglow emission, see lLazzati d2006h . 

GRBs and their afterglows can also be used to place 
constraints on Quantum Gravity (QG) because of a birefrin- 
gence effect on photon propagation, caused by the difference 
in light velocity for the two states of circular polarisation 
( Gambini & Pullinlll999|). Limits have b een obtained using the 
UV/optical afterglows dFan et al.ll2007t) . However, a definitive 
detection of the polarisation of the y-ray prompt emission will 
provide a much better constraint on models of QG. 



8. Conclusions 

The Spectrometer aboard INTEGRAL, SPI, has been used 
to measure the level of polarisation of the intense burst 



GRB 041219a which was detected by IBAS. The predicted in- 
strument response was obtained by Monte-Carlo simulations 
using the GEANT 4 mass model. Our results over several en- 
ergy ranges and two time intervals are consistent with a po- 
larisation signal of 60 + 35% which is a low level of signifi- 
cance (~ 2a-). The level of polarisation was calculated to be 
- 63;^3g% at an angle 70;^[j degrees for the 66 second 
time interval in the energy range 100-350 keV. The degree of 
polarisation was also constrained in the brightest 12 s of the 



GRB and a value of 96^-^q% at an angle of 60 degrees over 
the same energy range was obtained. Despite extensive analy- 
sis and simulations, we could not exclude a systematic effect 
that could mimic the weak polarisation signal. The polarisation 
fraction is within the range of the lower limits o btained from 



BATSE data for GRB 930131 and GRB 960924 dWilhs et al 
I2OO5,) . and also the value o f 41+^ ^ % obtained by RHESSI for 
GRB 021206 dWigger et alJl2004l) . 



As reviewed in ^ there are a number of model predictions 
available to explain the GRB observations. A significant degree 
of polarisation can be produced in GRBs by either synchrotron 
emission or by inverse Compton scattering. The level of po- 
larisation produced by synchrotron emission can be as high as 
70%. For Compton Drag, the condition FOj < 5 must be sat- 
isfied. In the case of GRB 041219a, F^^ ~ 3.3 and hence this 
process can explain significant y-ray polarisation. In some in- 
terpretations of GRB spectra, there can be a contribution from 
both the Compton and synchrotron processes. 
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